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Some features of hot and dense gas of quarks which are considered as the quasi-particles of the model Hamiltonian 
with four-fermion interaction are studied. Being adapted to the Nambu-Jona-Lasinio model this approach 
allows us to accommodate a phase transition similar to the nuclear liquid-gas one at the proper scale and to 
argue an existence of the mixed phase of vacuum and normal baryonic matter as a plausible scenario of crural 
symmetry (partial) restoration. Analyzing the transition layer between two phases we estimate the surface 
tension coefficient and discuss the possibility of quark droplet formation. 

-C 

Oh' Understanding in full and describing dependably the critical phenomena (chiral and deconfinement phase 

transitions) in QCD is still elusive because of a necessity to have the corresponding efficient non-perturbative 
methods for strongly coupled regime. For the time being such studies are pursued by invoking diverse effective 
tIh ■ models. The Nambu-Jona-Lasinio(NJL)-type models are certainly playing the most advanced role in this 

analysis pQ. This approach deals with the multi-fermion interactions in lieu of a gluon field QCD dynamics 
and does not incorporate the property of confinement. At the same time it is quite successful in realizing the 
^ ' spontaneous breakdown of chiral symmetry and its restoration at nonzero temperatures or quark densities. 

These and some related items are discussed in this paper inspired by well known and fruitful idea about 
the specific role of surface degrees of freedom in the finite fermi-liquid systems and in considerable extent by 
our previous works [2] and [3] in which the quarks were treated as the quasi-particles of the model Hamiltonian 
and the problem of filling up the Fermi sphere was studied in detail. Under such a treatment an unexpected 
singularity (discontinuity) of the mean energy functional as a function of the current quark mass was found. 
In the particular case of the NJL model the existence of new solution branches of the equation for dynamical 
quark mass as a function of chemical potential have been demonstrated and the appearance of state filled up 
with quarks which is almost degenerate with the vacuum state both in the quasi-particle chemical potential and 
in the ensemble pressure has been discovered. 

' Here we are going to study the quark ensemble features at finite temperature and fixed baryonic chemical 

potential and to analyse the first order phase transition which takes place in such a system of free quasi- 
particles. Analysis is performed within the framework of two approaches which are supplementary, in a sense, 
albeit giving the identical results. One of those approaches, based on the Bogolyubov transformations, is 
especially informative to study the process of filling the Fermi sphere up because at this point the density of 
quark ensemble develops a continuous dependence on the Fermi momentum. It allows us to reveal an additional 
structure in the solution of gap equation for dynamical quark mass just in the proper interval of parameters 
characteristic for phase transition and to trace its evolution. The result is that a quark ensemble might be 
found in two aggregate states, gas and liquid, and the chiral condensate is partially restored in a liquid phase. 
In order to make these conclusions easily perceptible we deal with the simplest version of the NJL model (with 
one flavor and one of the standard parameter sets) and, actually, do not aim to adjust the result obtained with 
well-known nuclear liquid-gas phase transition. Besides, it seems our approach might be treated as a sort of 
microscopic ground of the conventional bag model and those states filled up with quarks are conceivable as a 
natural 'construction material' for baryons. 

Now as an input to start with we remind the key elements of approach developed. The corresponding 
Hamiltonian includes the interaction term taken in the form of a product of two coloured currents located in 
the spatial points x and y which are connected by a form-factor and its density reads as 

H = -g(i 7 V + im)q - qt a ^q f dyq't b ^q' (A a A'»), (1) 
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where q = q(x), q = g(x), q' = q(y), q 1 = q(y) are the quark and anti-quark operators, 

q a i (x) = J (21^1)1/2 [a(p,s,c)u ai (p,s,c)e lpx + 6+(p, s, c)u m (p, s, c)e~ lpx ] , (2) 

p 2 = p 2 — m 2 , is the colour index, a is the spinor index in the coordinate space, a + , a and & + , 6 are the 
creation and annihilation operators of quarks and anti-quarks, a|0) = 0, b\0) = 0, |0) is the vacuum state of free 
Hamiltonian and m is a current quark mass. The summation over indices s and c is meant everywhere, the index 
s describes two spin polarizations of quark and the index c plays the similar role for a colour. t a = A a /2 are 
the generators of SU(N C ) colour gauge group, the Hamiltonian density is considered in the Euclidean space and 
7 M denote the Hermitian Dirac matrices, fx, v = 1,2,3,4. (A®A'„) stands for the form-factor of the following 

form (A^A'£) = S ab j^zz^ [/(x — y)S ll „ — J fiu (x — y)], where the second term is spanned by the relative distance 
vector and the gluon field primed denotes that in the spatial point y. The effective Hamiltonian density ([TJ 
results from averaging the ensemble of quarks influenced by intensive stochastic gluon field A^, see Ref. [2]. 
For the sake of simplicity in what follows we neglect the contribution of the second term of Eq.([l]). The ground 
state of the system is searched as the Bogolyubov trial function composed of the quark-anti-quark pairs with 
opposite momenta and with vacuum quantum numbers, i.e. 

\a) = T |0> , T = 11,,, exp{<^ [a+(p, s)b+(-p, s) + o(p, s)b(-p, s)}}. (3) 

In this formula and below, in order to simplify the notations we refer to one compound index only which 
means both the spin and colour polarizations. The parameter tp(p) which describes the pairing strength is 
determined by the minimum of mean energy E = (a\H\a). By introducing the 'dressing transformation' we 
define the creation and annihilation operators of quasi-particles as A — TaT" 1 , B + = Tb + T~ x and for fermions 
= T' . Then the quark field operators are presented as 

q{x) = J (2^3 (2| P !|)V 2 [ ^ (P ' S) U <P'4 elPX + B+ (P' s ) ^(P' s ) e ^ PX ] • 
9_(X)= / ( 2 t P )3 (2| P !|)V 2 [^ + (P^)^(P^)^ PX + 5(P^)^(P^)^ PX ] , 

and the transformed spinors U and V are given by the following forms U(p, s) = cos(^)tt(p, s) — s'm((p)v(— p, s), 
V(p,s) = sin(<ys)w(— p, s) + cos(<p)w(p, s) where U(p, s) = U + (p, 5)74, V(p,s) — V^ + (p, 5)74 are the Dirac 
conjugated spinors. 

In Ref. [3] the process of filling in the Fermi sphere with the quasi-particles of quarks was studied by 
constructing the state of the Sletter determinant type \N) = Il|p|<p F s ^ + (-f > ! S)\ a ) which possesses the minimal 
mean energy over the state | TV) . The polarization indices run through all permissible values here and the quark 
momenta are bounded by the limiting Fermi momentum Pp ■ The momenta and polarizations of states forming 
the quasi-particle gas are marked by the capital letters similar to above formula and the small letters are used 
in all other cases. 

As it is known the ensemble state at finite temperature T is described by the equilibrium statistical operator 
p. Here we use the Bogolyubov-Hartree-Fock approximation in which the corresponding statistical operator is 
presented by the following form 

P = 5 , Z a = Tr {e-^ H -^} , (4) 

where an approximating effective Hamiltonian H app is quadratic in the creation and annihilation operators of 
quark and anti-quark quasi-particles A + , A, B + , B and is defined in the corresponding Fock space with the 
vacuum state \a) and f3 — T^ 1 . There is no need to know the exact form of this operator henceforth because all 
the quantities of our interest in the Bogolyubov-Hartree-Fock approximation are expressed by the corresponding 
averages (a density matrix) n(P) = Tr{p A+(P; S)A(P; S)}, n(Q) = Tr{p S+(Q; T)B{Q; T)}, which are found 
by solving the following variational problem. The statistical operator p is determined in such a form in order 
to have at the fixed mean charge 

g 4 - Tr{p Q 4 } = V 2N C J -0^ [n(p) - n(p)] , (5) 

where Q4 = J [A + (p)A(p) + B(p)B + (p)] (for the diagonal component of our interest here, Q4 = 

— J dxqij4q) and fixed mean entropy 

S = —Trip In p\ = 

= -V2N C Hp) lnn(p) + (1 - n(p)) ln(l - n(p)) + n(p) lnft(p) + (1 - n(p)) ln(l - n(p))} , 
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(S = — hip) the minimal value of mean energy of quark ensemble E = Tr{pH}. The definition of mean charge 
is given here up to the unessential (infinite) constant coming from permuting the operators BB + in the charge 
operator Q4. It is reasonable to remind that the mean charge should be treated in some statistical sense because 
it characterizes quark ensemble density and has no colour indices. 

Calculating the corresponding matrix elements leads to the following result for mean energy density per one 
quark degree of freedom (the details can be found in [4]) w = £/2N c , £ = E/V where E is a total ensemble 
energy 



dp . f dp „ r 

-\PM + / ,r, \pi\cos9[n{p) + n(p) - 1\ 



(27r) 3 '^' J (2tt) 3 



(7) 



(up to the constant unessential for our consideration here). Here the following denotes are used p = |p|, q = |q|, 
9 — 2ip, 8' — 6(q), I — I(p + q) and the angle 9 m (p) is determined by sin6* m = m/|p 4 |. It was quite practical 
to single out the colour factor in the four-fermion coupling constant as G = 2G/N C . It is of important to notice 
that the existence of such an angle stipulates the discontinuity of mean energy functional mentioned above and 
found out in [2] 

We are interested in minimizing the following functional Q = E — fi Q4 — T S where p and T are the 
Lagrange factors for the chemical potential and temperature respectively. The approximating Hamiltonian 
^app is constructed simply by using the information on E — p Q4 of presented functional (see, also below). For 
the specific contribution per one quark degree of freedom / — F/2N C , F = il/V we receive 



/ = 



/ (2^3 cos6 ( n + Kn-n)]+ J J^TaM ~gJ sin (9 - 9 m ) (n + n - 1) x 

d f d ^ 

— %r sin (6 1 - 9') in' +fi'-l)I + T / [n\nn+ (1 - n) ln(l - n) + nhxn + (1 - fl) ln(l - ti)] . 

(27r) d J (2ny 

Here the primed variables correspond to the momentum q. The optimal values of parameters are determined 
by solving the following system of equations (df/d6 — 0, df/dn = 0, df/dn = 0) 

\p 4 \ sin 9 - M cos (9 - 9 m ) = , 

(9) 

\p±\ cos9-n + M sin (9 - 6 m ) - T In (ti^ 1 - l) = , 
\p 4 \ cos 9 + fj, + M sin (9 -6 m )-T In (ra^ 1 - l) = 

where we denoted the induced quark mass as 

M(p) =2G fj0z (1 - n' - n') sin (9' - 9'J I(p + q). (10) 

Turning to the presentation of obtained results in the form customary for mean field approximation we 
introduce a dynamical quark mass M q parameterized as sin (0 — 9 m ) — pfe, \P^\ — (p 2 + M^p)) 1 / 2 and 
ascertain the interrelation between induced and dynamical quark masses. From the first equation of system (|9|) 
we fix the pairing angle as sin# = pM / (\p4\\P^\) and making use of the identity 



(|p 4 | 2 - M m) 2 + M 2 p 2 = [p 2 + (M - m) 2 ] \p 4 
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find out that cos 9 = ± ^i^nr^r^ ■ F° r clarity we choose the upper sign 'phis'. Then, as an analysis of the NJL 
model teaches, the branch of equation solution for negative dynamical quark mass is the most stable one. Let 
us remember here that we are dealing with the Euclidean metrics (though it is not a principal point) and a 
quark mass appears in the corresponding expressions as an imaginary quantity. Now substituting the calculated 
expressions for the pairing angle into the trigonometrical factor expression sin (6 — 9 m ) = sm9 j^j — cos 9 t^t and 
performing some algebraic transformations of both parts of equation we define M q (p) — M(p) — m. And then 
the equation for dynamical quark mass (j!0[) is getting the form characteristic for the mean field approximation 



M = 2G /(03( 1 -^ / -^)^|^p + q) 



The second and third equations of system ([§]) allow us to find for the equilibrium densities of quarks and 
anti-quarks as n — [e& d P4 l~^) + n — [e 13 d P4 l+'') + and, hence, the thermodynamic properties of 

our system as well and, in particular, the pressure of quark ensemble P = —dE/dV. By definition we should 
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Figure 1. The chemical potential p (MeV) as a 
function of charge density Q\ = Qa/(ZV) (in the 
units of charge/fm 3 ). The factor 3 relates the den- 
sities of quark and baryon matter. The top curve 
corresponds to the zero temperature. The curves 
following down correspond to the temperature val- 
ues T = 10 MeV, ... , T = 50 MeV with spacing 
T = 10 MeV. 
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Figure 2. The ensemble pressure P (MeV/fm 3 ) 
as a function of charge density Q4 at temperatures 
T = MeV, ... , T = 50 MeV with spacing T = 10 
MeV. The lowest curve corresponds to zero tem- 
perature. The dashed curve shows the boundary 
of phase transition liquid-gas, see the text. 



calculate this derivative at constant mean entropy dS/dV = 0. This condition allows us, for example, to 
calculate the derivative d/i/dV. However, this way is not reliable because then the mean charge Q4 might 
change, and it is more practical to introduce two independent chemical potentials — for quarks /x and for 
anti-quarks \x (following formula for n with an opposite sign). In fact, it is the only possibility to obey both 
conditions simultaneously. It leads to the following definitions of corresponding densities n = [e@ d P4 l~^) + 
n = [e@ (1-^1+^) + On this way of description we are able even to treat some non-equilibrium states of 

quark ensemble (albeit with losing a covariance similar to the situation which takes place in electrodynamics 
while one deals with electron-positron gas). But here we are interested in the particular case of p, = p. Then 
the corresponding derivative of specific energy dw/dV might be presented as 



dw 
dV 



dp 



dn dfi 
djldV 



dn dp 
dfidV 



|p 4 |cos0-2Gsin(0-0 m ) 



dq 

(2^)3 



1)/ 



(11) 



Now expressing the trigonometric factors via dynamical quark mass and exploiting Eq. (j!0[) we obtain the ensem- 
ble pressure as P = -| - V2N C J + \Pi\- The requirement for mean charge conservation 

= 9± + V2N C J ^3 f - Ifffy = P rovides us with an equation which interrelates the deriva- 

tives dp/dV, dp/dV. Apparently, the regularized expressions for mean charge of quarks and anti-quarks are 
meant here. Dealing in a similar way with the requirement of mean entropy conservation, dS/dV — 0, we 



receive another equation as 



Tln(r 



1) 



"A 4 ' 



s _ _ r dp 

2N C V 2 J (2tt)3 

P A \ and Tlnih- 1 - 1) = 



ln(n 



1 \ dndf^, f 
L ) dudV T J 



dp 



ln(f 



Substituting here 



d^dV ' J (2tt) 3 

IP4I after simple calculations (keeping in mind that ji = fi 



and the charge conservation) we have that J jj^jj \djldv 



dp ( dn dn , dn dn \ i p i 

— ' — — + dp.-dV J - ' 
Q±l± 



ST 
2N C V 2 



2N C V 2 ■ 



Finally it leads 



for the pressure to the following expression P = — y + ^-T. _|_ !d±J± ( f course, the thermodynamic potential is 
= —PV). At small temperatures the anti-quark contribution is negligible and thermodynamic description 
can be grounded on utilizing one chemical potential \x only. If the anti-quark contribution is getting intrinsic 
the thermodynamic picture becomes complicated due to the presence of chemical potential fl with the condition 
jl — fi imposed. In particular, at zero temperature the anti-quark contribution is absent and we might receive 
P = —£ -\- /1 Pq where /i = [Pp + M^(Pp)] 1 / 2 , Pp is the Fermi momentum and p q = N/V is the quark ensemble 
density. 

For clarity, we consider the NJL model in this paper, i.e. the correlation function behaves as the <5-function 
in coordinate space. It is a well known fact that in order to have an intelligent result in this model one needs to 
use a regularization cutting of momentum integration in Eq. ([S]). We adjust the standard set of parameters [5] 
here with |p| < A, A = 631 MeV, m = 5.5 MeV and GA 2 /(2tt 2 ) = 1.3. This set of parameters at n = 0, n = 0, 
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Figure 3. The fragments of isotherms in Fig. [IJ 
see text. Chemical potential fj, (MeV) as a func- 
tion of pressure P (MeV/fm 3 ). The top curve cor- 
responds to the zero isotherm and following down 
with spacing 10 MeV till the isotherm 50 MeV (the 
lowest curve). 
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Figure 4. The surface tension coefficient u s in 
MeV as a function of parameter c (C = c rf) for the 
curve of stable kinks (with r\ < 1.2). 



T = gives for the dynamical quark mass M q = 335 MeV. Besides, it may be shown that the following form of 
ensemble energy is valid at the extremals of functional © 

E = E vac + 2N C V J j^-\P 4 \(n + n) , E vac = 2N c vJ J^Qp*] - \P*\) + 2N - V ^ • 

It is easy to understand that this expression with the vacuum contribution subtracted looks like an energy of a 
gas of relativistic particles and anti-particles with the mass M q and coincides identically with that calculated 
in the mean field approximation. 

Thus, we determine the density of quark n and anti-quark n quasi- particles at given parameters it and T 
from the second and third equations of system (|9]). From the first equation we receive the angle of quark and 
anti-quark pairing 6 as a function of dynamical quark mass M q which is handled as a parameter. The evolution 
of chemical potential as a function of charge density Q 4 = Qi/(3V) (in the units of charge//™ 3 ) with the 
temperature increasing is depicted in Fig. [T] (factor 3 connects the quark and baryon matter densities). The 
top curve corresponds to the zero temperature. The other curves following down have been calculated for the 
temperatures T = 10 MeV, ... , T = 50 MeV with spacing T = 10 MeV. As it was found in Ref. [3] the chemical 
potential at zero temperature increases first with the charge density increasing, reaches its maximal value, then 
decreases and at the densities of order of normal nuclear matter densit}0, p q ~ 0.1 6/ /to 3 , becomes almost equal 
to its vacuum value. Such a behaviour of chemical potential results from the fast decrease of dynamical quark 
mass with the Fermi momentum increasing. It is clear from Fig. Q] that the charge density is still a multivalued 
function of chemical potential at the temperature slightly below 50 MeV. Fig. [2] shows the ensemble pressure 
P (MeV/fm 3 ) as the function of charge density Q4 at several temperatures. The lowest curve corresponds 
to the zero temperature. The other curves following up correspond to the temperatures T = 10 MeV, ... , 
T = 50 MeV (the top curve) with spacing T = 10 MeV. It is interesting to remember now that in Ref. [5] the 
vacuum pressure estimate for the NJL model was received as 40 — 50 MeV/fm 3 which is entirely compatible 
with the results of the conventional bag model. Besides, some hints at instability presence (rooted in the 
anomalous behavior of pressure dP/dn < 0) in an interval of Fermi momenta has been found. Fig. [3] shows the 
fragments of isotherms of Fig. Q] [2] but in the different coordinates (chemical potential — ensemble pressure). 
The top curve is calculated at the zero temperature, the other isotherms following down correspond to the 
temperatures increasing with spacing 10 MeV. The lowest curve is calculated at the temperature 50 MeV. The 
Fig. [3] obviously demonstrates a presence of the states on isotherm which are thermodynamically equilibrated 
and have equal pressure and chemical potential (see the characteristic Van der Waals triangle with the crossing 
curves). The calculated equilibrium points are shown in Fig. [2] by the dashed curve. The intersection points 
of dashed curve with an isotherm are fixing the boundary of gas — liquid phase transition. The corresponding 
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straight line P — const which obeys the Maxwell rule separates the non-equilibrium and unstable fragments of 
isotherm and describes a mixed phase and appropriate critical temperature for the parameter we are using in 
this paper turns out to be T c ~ 45.7 MeV with the critical charge density as Q4 ~ 0.12 charge/fm 3 . Usually 
the thermodynamic description is grounded on the mean energy functional which is the homogeneous function 
of particle number like E — N f(S/N,V/N) (without vacuum contribution). It is clear that such a description 
requires the corresponding subtractions to be introduced, however, this operation does not change the final 
results considerably. It was argued in Ref. i3j that the states filled up with quarks and separated from the 
instability region look like 'natural construction material' to compose the baryons and to understand the existing 
fact of equilibrium between vacuum and octet of stable (in strong interaction) baryonfl 

Apparently, our study of the quark ensemble thermodynamics produces quite reasonable arguments to 
propound the hypothesis that the phase transition of chiral symmetry (partial) restoration has already been 
realized as the mixed phase of physical vacuum and baryonic matteo However, it is clear our quantitative 
estimates should not be taken as ones to be compared with, for example, the critical temperature of nuclear 
matter which has been experimentally measured and equals to 15 - 20 MeV. Besides, the gas component (at 
T = 0) has nonzero density (as 0.01 of the normal nuclear density) but in reality this branch should correspond 
to the physical vacuum, i.e. zero baryonic densitjQ. In principle, an idea of global equilibrium of gas and liquid 
phases makes it possible to formulate the adequate boundary conditions at describing the transitional layer 
arising between the vacuum and filled state and to calculate the surface tension effects. 

The idea advanced would obtain substantial confirmation if it becomes possible to claim an evidence of 
existing the transition layer at which the ensemble transformation from one aggregate state to another takes 
place. As it was argued above the practical parameter for describing an uniform phase (at a given temperature) 
is the mean charge (density) of ensemble. Thus, one can reconstruct all other characteristics, for example, a 
chiral condensate, dynamical quark mass, etc. Analyzing the transition layer at zero temperature we assume 
the parameters in the gas phase are approximately the same as at zero charge p g = 0, i.e. as in the vacuum 
(ignoring the negligible distinctions in the pressure, chemical potential and quark condensate). Then dynamical 
quark mass obtained has maximal value and for the parameter choice of the NJL model it is M = 335 MeV. 
From the Van der Waals diagram one may draw out that the second (liquid) phase being in equilibrium with 
the gas phase develops the density pi — 3 x 0.185 ch/fm 3 . The detached factor 3 here relates the magnitudes 

of quark and baryon densities. The quark mass in this phase is approximately M~ 70 MeV (we are dealing 
further with the simple one-dimensional picture). 

Usually an adequate description of heterogeneous states can be developed basing on the mean field approxi- 
mation [7j, specifically for our case, by dealing with the corresponding effective quark-meson Lagrangian (a sort 
of the Ginzburg-Landau functional) 

I ! 7T1 2 

£ = -q(d + M) q - - (d^a) 2 - U(a) - - F^F^ - - g a qq a + ig v q 7 „ q , (12) 

2 

where F^ = d^V v — <9„V^, U (a) — ^-er 2 + |cr 3 + f er 4 , a is the scalar field, is the field of vector mesons, m CT , 
m v are the masses of scalar and vector mesons and g ai g v are the coupling constants of quark- meson interaction. 
The U(a) potential includes the nonlinear terms of sigma-field interactions up to the fourth order, for example. 
For the sake of simplicity we do not include the contribution coming from the pseudoscalar and axial-vector 
mesons. 

We are not going beyond well elaborated (and quite reliable) one loop approximation (fT2]) , although recently 
the considerable progress was reached in scrutinizing the non uniform quark condensates by utilizing the powerful 
methods of exact integration [8]. We believe it is more practical to adjust phenomenologically the parameters 
of effective Lagrangian being guided also by transparent physical picture. It is easy to see that handling one 
loop approximation actually we have the Walecka model [9] but applied for the quarks. In what follows we are 
working with the notations of that model hoping it does not lead to the misunderstandings. In the context 
of our deliberation Eq. (TT21 can be interpreted in the following way. Each phase, in some extent, might be 
considered as an excited state as to its relation with another phase which requires an additional (besides a charge 
density) set of parameters just as the meson fields for describing and those fields characterize the measure of 
deviation from the equilibrium state. Then the key question becomes whether it is possible to adjust the effective 
Lagrangian parameters of (|12jl in order to obtain the solutions in which the quark field interpolates between 
the quasi-particles in the gas (vacuum) phase and in the quasi-particles of the filled up state. The density of 
ensemble of the filled up states should asymptotically approach an equilibrium value of p\ and turn to zero value 
in the gas phase (vacuum). 

2 The chiral quark condensate for the filled up state discussed develops the quantity about (100 MeV) 3 (at T = 0), see [3J, that 
demonstrates the obvious tendency of restoring a chiral symmetry. 

3 Indireet confirmation of this hypothesis one could see, for example, in the existing degeneracy of excited baryon states Ref. [6]. 

4 Similar uncertainty is present in the other predictions of chiral symmetry restoration scenarios, for example, it stretches from 
2 to 6 units of normal nuclear density. 
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Taking the parameterization of the potential U{&) as b a = 1.5 m 2 (g a /M), c a = 0.5 m 2 (g a /M) 2 we come 
to the sigma model and the choice b = 0, c = results in the Walecka model. As to the application for nuclear 
matter the parameters b and c demonstrate essentially the model dependent character and arc different from 
the parameter values of the sigma model. They are phenomenologically adjusted with requiring an accurate 
description of the saturation property. On the contrary, for the quark Lagrangian (|12[) we could intuitively 
anticipate some resemblance with the sigma model and, hence, introduce two dimensionless parameters 77, C 
as b — r] b a , c — ( 2 c a which characterize some fluctuations of the effective potential. Then the scalar held 



potential is presented by the following form U(a) = ^f-jf? {^^§ T 
fields are defined by the following system of the stationary equations 



4^-770- + £ z a< 



ml a 



A a 
AV-mlV 



b a z +c er 
: ~9v P , 
(V+ M) q = (E — g v V) q 



9a Ps 



The meson and quark 



(13) 



where M= M +g a o~ is the running value of dynamical quark mass, E stands for the quark energy and V = —1V4,. 
The density matrix describing the quark ensemble at T = has the form £(x) — J Ff j^3q P (x)q p (x) where p 
is the quasi-particle momentum and the Fermi momentum Pp is defined by the ensemble chemical potential. 
The densities p s , p in the right hand sides of equations (|13|) equal (by definition) to p s (x) = Tr{£(x),l}, 
p(x) = Tr {£(x), 74}. Here we confine ourselves to the Thomas-Fermi approximation while describing the 
quark ensemble. Then the densities in which we are interested in are given with some local Fermi momentum 

P F (x) as p = 7 J Pf = g^F| p s = 7 J Pf where 7 is the quark gamma factor (7 = 2N c N f , N c 



is the number of colours, Nf number of flavours), E = (p 2 + M ) 1 ^ 2 - By definition the ensemble chemical 
potential does not change and it leads to the situation in which the local value of Fermi momentum is defined 

* 2 

by the running value of dynamical quark mass and vector field as p. = M = g v V + (P F + M ) ' . The 
details of tuning the Lagrangian parameters (TT2t can be found in [10] . The point of our attraction here is the 

surface tension coefficient [IT] u s = Air r 2 , dx £(x) — ^ p(x) 



here £1 is the energy density in the liquid 

phase. The parameter r Q is discussed below. In the Thomas-Fermi approximation £{x) — 7 J Pf ^ [p 2 + M 

(a;)] 1 / 2 + ^g v p(x)V(x) — ^g a p s (x)a(x). The surface tension coefficient u s in MeV for the curve of stable kinks 
(see the details in Ref. [TU]) with parameter 77 < 1.2 as the function of another parameter c (£ = c 77) is depicted 
in Fig. SI 

Above results lead us to put the challenging question about the properties of finite quark systems or droplets 
of quark liquid which are in equilibrium with the vacuum state. As a droplet here we imply the spherically- 
symmetric solution of the equation system ()13[) for o~(r) and V(r) with the obvious boundary conditions cr'(0) = 
and V'(0) = in the origin (the primed variables denote the first derivatives over r) and rapidly decreasing 
at the large distances a — >■ 0, V — » when r — >■ 00. Fig. [5] shows the set of solutions (cr-field (MeV)) of the 
equation system (fi"3j) at number of flavors Nj = 1. Fig. [6] presents the corresponding distributions of ensemble 
density p (ch/fm 3 ). The Table 1 exhibits the results of fitting the density p(r) with the Fermi distribution 

Pf(t) = 1+e (ff- r )/b where po is the density at the origin, Rq is the mean size of the droplet and the parameter 

b determines the thickness of surface layer t — 41n(3)6. Besides, the coefficient ro which is included in the 

1/3 

definition of surface tension coefficient, Rq = r$N q is also presented together with the characteristic values of 
the cr-meson mass and the coefficient 77 at which all this values were obtained. 

Table 1 . Results of fitting by the Fermi distribution (Nf = 1). 





po (ch/fm 3 ) 


Ro (H 


b (fin" 1 ) 


t (fm) 


r (fm) 


m a (MeV) 


V 


15 


0.34 


1.84 


0.51 


2.24 


0.74 


351 


0.65 


43 


0.43 


2.19 


0.52 


2.28 


0.75 


384 


0.73 


159 


0.46 


4.19 


0.52 


2.29 


0.77 


409 


0.78 


303 


0.47 


5.23 


0.52 


2.29 


0.78 


417 


0.795 


529 


0.47 


6.37 


0.52 


2.27 


0.79 


423 


0.805 


742 


0.47 


7.15 


0.52 


2.27 


0.79 


426 


0.81 



The curves plotted in the Fig. [5] and results of Table 1 allows us to conclude that the density distributions at 
N q > 50 correspond quite well to the data typical for the nuclear medium. The thicknesses of transition layers 
are also similar. The coefficient ro with the factor 3 1 / 3 included is in full correspondence with nuclear one. The 
values of the a-meson mass turn out to be quite reasonable as well. Although at small quark numbers in the 
droplet the corresponding behaviors become essentially different. We know experimentally that in the nuclear 
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Figure 5. a- field (MeV) as a function of the 
distance r (fm) for several solutions of the equa- 
tion system fj 13() which are characterized by the 
net quark number N q written to the left of each 
curve. 



Figure 6. Distribution of the quark density p 
(ch/fm 3 ) for the corresponding solutions presented 
in the Fig. El 



matter one can observe some increase of the ensemble density which is quite considerable for the Helium and is 
much higher than the normal nuclear matter density for the Hydrogen. One may criticize us in this point because 
working within the Thomas-Fermi approximations becomes hardly justified at the small number of quarks and 
it is necessary to handle the solution of equation system ((TBI . However, fortunately, the exploration we are 
interested in has been performed in the chiral soliton model of nucleon [12] . It has been demonstrated there that 
adding the contributions of pseudo-scalar and axial- vector helds to the Lagrangian (fT2|) leads to reasonably good 
description of nucleon and A. The interesting remark here is that the soliton solutions obtained in [12] could be 
interpreted as a confluence of two kinks. Each of those kinks develops the restoration of chiral symmetry in a 
sense that the scalar field is approaching its zero value at the distance ~ 0.5 fm from the kink center. Actually, 
one branch corresponds to the solution with the positive value of the dynamical quark mass and another branch 
presents the solution with negative dynamical quark mass (in three-dimensional picture the pseudo-scalar fields 
appears just as a phase of chiral rotation from positive to negative value of quark mass). 

Acknowledgements. We are grateful to the meeting organizers and personally Professor V. Skalozub for a 
hospitality and excellent conditions for efficient work. 
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